Static wetting on deformable substrates, from liquids to soft solids 
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Young's law fails on soft solid and liquid substrates where there are substantial deformations near 
the contact line. On liquid substrates, this is captured by Neumann's classic analysis, which provides 
a geometrical construction for minimising the interfacial free energy. On soft solids, the total free 
energy includes an additional contribution from elasticity. A linear-elastic model incorporating an 
out-of-plane restoring force due to solid surface tension was recently shown to accurately predict the 
equilibrium shape of a thin elastic film due to a large sessile droplet. Here, we extend this model 
to find substrate deformations due to droplets of arbitrary size. While the macroscopic contact 
angle matches Young's law for large droplets, it matches Neumann's prediction for small droplets. 
The cross-over droplet size is roughly given by the ratio of the solid's surface tension and elastic 
modulus. At this cross-over, the macroscopic contact angle increases, indicating that the substrate 
is effectively less wetting. For droplets of all sizes, the microscopic behaviour near the contact line 
follows the Neumann construction giving local force balance. 
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I. INTRODUCTION 

Wetting is a fundamental physical process with far- 
ranging applications [1-4] . In the absence of long-range 
interfacial forces, our understanding of wetting typically 
centres upon two key results: Young's law for wetting 
on rigid substrates and Neumann's triangle for wetting 
on liquid substrates [1, 5-7]. These are shown schemat- 
ically in Figure 1. On a rigid substrate, the total free 
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FIG. 1. Wetting on (a) rigid and (b) liquid substrates as 
described by Young and Neumann, respectively. 

energy is minimized when 7s„ = ')si+liv cos 9. Here, 7si, 
7su and 7;^, are the intensive free energies for the solid- 
liquid, solid-vapour and liquid-vapour interfaces, respec- 
tively [1]. This result is often visualized in terms of an 
in-plane balance of surface tensions, whose magnitudes 
are given by the interfacial energies, acting at the con- 
tact line. This interpretation of Young's law is not rig- 
orously correct since solid surfaces can have contibutions 
to the surface tension beyond the interfacial energy [8]. 
On a liquid substrate, the three interfaces arrange so that 
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surface tensions balance in- and out-of-plane at the con- 
tact line. This configuration, called Neumann's triangle, 
uniquely determines the angles between the phases at 
the contact line [6, 9]. For liquids, the equivalence of sur- 
face tensions and interfacial energies is not problematic. 
For small drops, where the effect of gravity is unimpor- 
tant, droplets on a liquid substrate then assume lentic- 
ular forms, such as that shown in Figure 1(b), where 
the interfaces between phases take the form of spherical 
caps[l]. 

Young's law and Neumann's triangle are useful in ex- 
plaining many wetting phenomena. However, in a sense 
they are two extremes in a continuum of wetting prob- 
lems: Young's law gives the behaviour in the limit of 
an infinitely hard substrate, while Neumann's triangle 
gives the result for an infinitely soft [i.e. fiuid) substrate. 
Wetting on soft solid substrates should then exhibit a 
spectrum of intermediate behaviour. This liquid-solid 
crossover is suggested by several recent results. Exper- 
iment, analytical theory and molecular-dynamics simu- 
lations, suggest that the shape of a soft substrate near 
a contact line does not depend on the substrate elastic 
modulus [10, 11]. Eslami and Elliott[12] showed that 
observations of condensation on soft substrates can be 
explained by treating soft surfaces as fiuids. Mora et al. 
[13] showed that the Rayleigh-Plateau instability can oc- 
cur in thin elastic strands. 

In this paper, we generalize Jerison et al.'s solution for 
a large droplet on a thin solid substrate to droplets of 
arbitrary radius. For large droplets, we recover Young's 
law for the macroscopic contact angle. For small droplets, 
the macroscopic behaviour is given by a lenticular shape, 
such as that seen in Figure 1(b), as predicted by Neu- 
mann's analysis for wetting on liquid substrates. The 
cross-over droplet size is approximately the ratio of the 
solid's surface tension and elastic modulus. For droplets 
of all sizes, the microscopic behaviour near the contact 



line follows a Neumann construction so as to give local 
force balance. 



II. DEFORMATION OF AN ELASTIC SLAB 
WITH SURFACE TENSION 

Recently, Jerison et al. [10] measured the deforma- 
tion of a thin film of soft silicone gel under a large water 
droplet. Their results showed that the surface profile was 
well-matched by a linear-elastic theory that included the 
out-of-plane restoring force due to the substrate's surface 
tension. The theory they presented was two-dimensional, 
and thus applicable to very large droplets. However in 
order to apply the theory to droplets of a finite radius, 
we need to extend their results to three dimensions. To 
this end, we solve the elastic equations in cylindrical co- 
ordinates for a uniform substrate of thickness h and of 
infinite horizontal extent. 



A. Purely elastic deformations 

We take a substrate which is pinned to a rigid sur- 
face at z = with displacements there being zero. In 
cylindrical-polar coordinates, the governing equations are 
the steady-state Navier equations 

(l-2^)fvV-^)+|-(V.u) = 
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where we have assumed cylindrical symmetry, Ur and u^ 
are the r and z components of the displacement u, and 
v is the Poisson's ratio[14]. The stress is related to the 
displacements by 
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where E is Young's modulus. Because of the radial 
symmetry, we take Hankel transforms of the stress 
and displacement fields, giving Ur{s,z) = tli[ur{r, z)], 
Uz{s,z) = Ho[uz(r, 2;)], arz{s,z) = Hi [cTrz (r, z)] and 
azz{s,z) = Hol'Tzzii'T z)]. Here Hq and Hi are Hankel 
transforms of order and 1 respectively [15]. By follow- 
ing the same solution method used Xu et al. [16], we 
then obtain 



ai{s,h) = Qij{s,h,z)uj{s,z), 



(3) 



where ai{s,z) — {arz{s, z),azz{s, z)) and Ui{s,z) — 
(Ur{s,z),Uz{s, z)). The spring constant Qij{s,h,z) that 
relates surface tractions to substrate displacement is 
given by 
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where 

Pij{s,h) = 
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Knowing Qij enables the calculation of displacements 
throughout a substrate given forces applied to its sur- 
face, as it is directly related to the Green's function for 
the problem. For example, in the case of purely vertical 
forces acting on the surface, the surface displacements are 
given by Uz{r, h) = aQ^[Q~}{s, h, h)azz{s, h)], where 
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(7) 

As a useful check, Q~}{s, oo, oo) = 2(1 — i'')/{sE) which 
gives Terezawa's displacement solution for axisymmetric 
forces acting on a semi- infinite elastic substrate [15]. 



B. Elasticity with surface stresses 

The elastic response due to stresses applied at the free 
surface is given by Equation (3). However, to describe 
the deformation due to the wetting of a liquid droplet, we 
also need to include surface stresses - a generalization of 
surface tension for solids. The surface stress, T, is related 
to the surface energy 7 by the Shuttleworth equation 
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where Sij is the Kronecker delta, and e^j is the surface 
strain [2, 8]. T is generally anisotropic, however it does 
simplify under certain situations. For instance at a fluid- 
fiuid interface -g^ = 0, so the surface stress is isotropic 
and equal to the surface energy. Both T and 7 can then 
be called the 'surface tension' without ambiguity. For 
many isotropic solid materials the surface stresses are 
also approximately isotropic [17-19], and T and 7 are 
typically of similar magnitudes [18]. However they are 
not necessarily equal [19, 20] and we have to distinguish 
carefully between them. Here we follow Shuttleworth [8] 
in referring to T as the 'surface tension' with the under- 
standing that it represents the actual tension force at the 
surface of the solid. 



For tractability, we assume isotropic surface tensions, 
and that the sohd-vapour surface tension T^^, and the 
sohd-hquid surface tension T,,; arc the same and given by 
Tg. Later we discuss expected changes in more general 
cases. 

The hnearised surface tension force is ctt = 
TTs-^^ ( ^ " Or ' ) 2, where z is the unit vector in the 
vertical direction [18]. We include this in the force bal- 
ance at the surface [10] to find that 



ai{s,h) = QSij{s,h,z)uj{s,z), 



where 



QSij{s, h, z) = Qi]{s, h, z) + 
In particular, we have 




T,s 



(9) 



(10) 



sE 



(3-4iy)sinh(2s/i)-2sh "*" E 



(11) 



QS^^{s, h, h) tells us how the surface of the substrate 
responds to imposed forcings of wavelength 0(l/s), and 
we can use this to gain some interesting insights into 
the physics of the problem. For small wavenumbcrs, 
QS-^is,h,h) -^ Q-^{s,h,h). That is, for long- 
wavelength surface perturbations, the response of the 
substrate is purely elastic. On the other hand, for large 
s, QS~^{s,h,h) — > l/{Tss'^). That is, for short wave- 
length perturbations to the surface, the force due to sur- 
face tension dominates the substrate response. Analysis 
of Equation (11) shows us that the cross-over length scale 
is Tg/E, provided the substrate thickness h > Tg/E. 
Perturbations of lengthscale A ;3> Tg/E arc damped elas- 
tically, while perturbations of length scale A ^ Tg/E 
are damped by surface tension [21]. For thin substrates 
h < Tg/E, the substrate is less comphant due to the 
presence of the rigid bottom boundary, so the elastic re- 
sponse to perturbations is stronger, and the crossover 
lengthscale is reduced. As we shall see, the length scale 
Tg/E repeatedly emerges in substrate micro-deformation 
problems as a controlling influence on the behavioural 
response of the substrate. Because it represents the bal- 
ance between elasticity and capillarity we refer to it here 
as the elasto-capillary length, though it should be noted 
that there are several other length scales that have also 
been given this title [22]. 

At this point it is worth correcting a small error in the 
analysis of Jerison et al. [10]. When including the sur- 
face tension of the substrate in their model, they used 

(Tt = 2 — "j^f ' ■ The correct linearised surface tension 
force does not have the factor of 1/2. This can be de- 
rived from the fact that for a curved surface ctx = Tg/Cn 
where /C is the surface curvature and n is the normal 
vector to the surface [81. The surface-tension force is the 



linearised version of this expression. Importantly this 
correction does not affect Jerison et aUs claim that that 
a linearised elastic model matches the observed surface 
profile. Rather, it changes the fitted ratio of the liquid 
and solid surface tensions. 



III. SOLUTION FOR A HEMISPHERICAL 
DROPLET 

Using the results above, we can calculate the exact 
deformation of the surface caused by a hemispherical 
droplet of radius R. Assuming that all length scales 
in the problem are much bigger than the typical range 
of intermolecular forces, then the surface tension of the 
droplet appears as a line force, and the traction imposed 
by the droplet on the surface of the substrate is given by 



T{r,h)=-fiS{r-R)i-PiH{R 



(12) 



where 7/ is the liquid/vapour surface tension. Pi is the 
Laplace pressure in the droplet, 6{x) is the Dirac delta 
function, H{x) is the Heaviside step function, and z is 
the unit vector in the z direction [23]. The first term 
corresponds to the out-of-plane surface tension force of 
the droplet at the contact line. Here we have assumed 
that the liquid- vapor interface is oriented normal to the 
substrate in accordance with Young's law. We will re- 
lax this constraint in Section IV. The second term corre- 
sponds to the Laplace pressure force. To ensure mechan- 
ical equilibrium, the total force exerted by the droplet 
on the substrate is zero, thus the distributed Laplace 
pressure inside the droplet must balance the localised 
force of the droplet surface tension at the contact line[l]. 
For our hemispherical droplet, this means P; = 2"/i/R. 
Taking the Hankcl transform Ho[T(r, ft,)], we find that 
Ozz{s,h) — jiRJq{sR) — 2jiJi{sR)/s, where Ji{z) is the 
ith order Bessel function of the first kind. Therefore, us- 
ing Equation (9) and taking the inverse Hankel transform 
we find the vertical surface deformation of the substrate 
that results from the presence of the droplet: 
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The horizontal displacements Ur{r,z) can be similarly 
calculated. 



A. Large droplets on a thin substrate: R^ h 

There are three asymptotic limits that are of particu- 
lar use. Firstly, when the droplet is very large compared 
to the substrate thickness, we expect that the peak pro- 
file should approach the two-dimensional solution given 
by Jerison et al. [10]. If we define e = h/R, nondimen- 
sionalise s by setting s = sh, and move into the frame 



of reference of the peak by setting x - 
we can expand the results for u^ and u 



(r - R)/h. Then 
from above in 
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powers of e, to find that at leading order the surface dis- 
placements are 
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These arc the same as the results given by Jcrison et al. 
[10], after the correction of a small error in their anal- 
ysis noted above. Thus we recover the two-dimensional 
solution. 

Figures 2(a) and (b) demonstrate how the vertical 
and horizontal surface displacements tend to the two- 
dimensional solution as the droplet diameter increases. 
For smaller droplets, the shape of the wetting ridge is 
asymmetric, as has been observed experimentally [24]. 
As the droplet radius increases the shape of the wet- 
ting ridge becomes more symmetric. For R/h > 0(10), 
the two-dimensional solution provides a good estimate to 
Uz{r,h) near the contact line. However a significant dif- 
ference persists between Ur{r, h) and the two-dimensional 
theory up to relatively large droplet radii. This can 
explain the deviation between experiments and theory 
found by Jerison et al. [10]. They measured droplets with 
R/h ~ 0(10) and found excellent agreement between the 
two-dimensional theory and experiments for vertical dis- 
placements. However experimentally-observed horizon- 
tal displacements were asymmetric and shifted towards 
positive displacements relative to theoretical predictions. 
Our results suggest that this is simply due to the finite 
droplet size. 



B. Small droplets: i? < ft, i? < Ts/E 

Consider a sessile droplet of radius R <C Tg/E on a 
substrate of thickness h ^ R. Since the wavelength of 
applied surface stresses can be no bigger than than the 
droplet radius, then the arguments from the end of Sec- 
tion H imply that surface tension dominates the response 
of the substrate. As elasticity effectively drops out of the 
problem at these small lengthscales, the system is the 
same as that of a droplet of liquid sitting at the interface 
between two other fluids. Thus, the solution must simply 
be that of Neumann's classic three-fiuid contact problem 
shown in Figure 1(b), with ^si and 7^^ replaced by Tg. 
Here, the droplet takes a lenticular shape formed by the 
union of two spherical caps, with contact angles given by 
Neumann's triangle [9]. 



C. Droplets on a semi-infinite substrate: R <^ h, 

Ts/E < h 

In the case of an infinitely-thick substrate, as sh — > cxd, 
we find that Equation (11) simplifies considerably. We 
nondimensionalise using the droplet radius as a length 
scale, so that f = r/R and s = sR to find that the 
surface displacement is 
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This equation shows two interesting details. Firstly, 
the height of the peak has a linear dependence on 'yi/E, 
as predicted by many theoretical works [25], but is also 
a weak function of RE/Tg. Asymptotically we find that 
for R > Ts/E, 
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RE 
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+ const. (17) 



Figure 2(b) shows the peak height as a function of droplet 
radius, along with this leading order asymptotic result 
without the constant term. Evidently, the asympotic ex- 
pression is a convenient upper bound on the surface dis- 
placement of a substrate caused by a drop of radius R. 
This is in contrast to Terezawa's solution for the case of a 
circular line force on a semi-infinite substrate, where the 
substrate strain diverges at the contact line[15], and the 
two-dimensional models of Jerison et oL[10] and Long 
et o/.[21, 26] that predict a divergent peak height for 
semi-infinite substrates. Our asymptotic expression re- 
sembles the logarithmic divergence of peak height with 
R found by White [27]. However, while White's expres- 
sion depends on the range of intermolecular forces, our 
result only depends on drop radius and the elastocapil- 
lary lengths Ts/E and ■'yi/E. 



IV. DISCUSSION 
A. Does Young's law hold on soft substrates? 

Above, we have imposed a particular macroscopic con- 
tact angle 9 for the droplet, namely the value given by 
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FIG. 2. Substrate deformation with droplet size (a) out-of- 
plane displacements, and (b) in-plane displacements. The 
drop is to the left of the wetting ridge, h = 50/im, E — 
3kPa, Ts = 0.06N/m and u = 1/2. For R = 0.2 mm (blue 
dashed-dotted curve) and 0.4 mm (red dashed curve), the 
profile is computed using Equation (13) and its Ur analogue. 
The large droplet limit is computed using Equations (14) and 
(15) (black curve), (c) The height of the wetting ridge on a 
semi-infinite substrate. The dashed line shows the asymptote 
for RE/Ts » 1. 



Young's law. However, some experiments have suggested, 
this may not be the equilibrium contact angle on a soft 
substrate [28]. Therefore, we seek to determine the range 
of validity of Young's law using the model derived above. 
We calculate the equilibrium shape of a sessile droplet 
on a soft substrate by minimizing its free energy, written 
as i^ = Fsurf + Fei , where the contribution from surface 
energies. 



Fsurf = llAv + ilsl - 'ysv)As 



(18) 



with Aiy being the area of the liquid-vapour interface, 
and Asi being the area of the substrate-liquid interface. 
Note that we use the surface energies 7si and 7su here 
rather than the surface stresses. The elastic contribution 



F,i 



T.udA. 



(19) 



Note that we ignore the contribution to F of the line 
tension r. This is only expected to be significant when 
the size of a droplet ~ t/ji [1]. Typically r ~ lO^^^J/ni 
[1], so line tension can generally be ignored for droplets 
larger than a nanometre. 

The contact angle for a droplet is then found by min- 
imising this energy subject to the constraint that the 
total droplet volume is constant. The elastic energy 
starts to impact the equilibrium contact angle when 
Fei '-^ Fsurf and we can predict when this will occur 
by noting that for our hemispherical droplet, 



TrR'jiUz{R,h) 
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ru^{r,h)dr. (20) 



The first term corresponds to the work done by the sur- 
face tension of the droplet in pulling up the wetting ridge, 
while the second term corresponds to the work done by 
the Laplace pressure of the droplet. Each of these two 
expressions scale like jiRu^^R, h), and so using Equation 
(17) gives an upper bound on the magnitude of the elas- 

tic energy: F,i < ^\og{RE/Ts). Fsurf - liR^ so 
we find that the ratio F^i/ Fsurf ^ {li/RE) \og{Ts/ RE). 
This is small provided that R ^ Ts/E,ji/E, in which 
case the deformation of the substrate will have little in- 
fluence on the energetics of the droplet. Then the contact 
angle will be the same as it is for a rigid substrate. On 
the other hand, as the droplet radius reduces towards 
uiax.[Ts/E,^i/E], the growing influence of the substrate 
deformation will start to manifest itself as a change in 
the macroscopic contact angle. 

The key result here is that Young's law is recovered 
for droplets that are much larger than the two elasto- 
capillary lengths Tg/E and ji/E, i.e. ^sv = Isi + 
7/1, cos 9. An instructive example of this is given in Fig- 
ure 3(a). In this case, R ^ T s/E,^i/E. Therefore, 
macroscopically the droplet conforms to Young's law, 
with 9 = 90° as we assume 7^/ = ^sv However if we 
zoom in on the contact line, we see that the local angle 
between the liquid-air and substrate liquid interfaces de- 
viates substantially from 9. This is clearly demonstrated 
in Figure 3(b). Note that this second figure is calculated 
using Equation (13) and the vertical length scale is non- 
dimensionalised with ^1/ E. 

In order to demonstrate how Young's law holds for 
large droplets, and also how the contact angle changes as 
R approaches Ts/E and 'ji/E, we calculate the change 
in the macroscopic contact angle, 9, for an incompress- 
ible, thin substrate. The details of this calculation are 
given in Appendix A. Figure 4(a) shows the macroscopic 
contact angle for a droplet of volume V = 4/37tRq, with 
E = 3kPa, h = 20/im and 7; = 0.07N/m. 7^/ = 7s.„ 
so the macroscopic contact angle given by Young's law 
is 90°. We plot the curve for 6 for T^ = 0.08, 0.12 
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FIG. 3. Substrate deformation for a hemispherical 1 mm radius droplet of water on an incompressible elastic substrate of 
thickness h = 0.5mm with T^ = 0.05N/m and E = 3kPa. (a) Macroscopic view, (b) Substrate displacement underneath 
the droplet. Here, the vertical axis is scaled by Ji/E. (c) Close up of the tip of the wetting ridge, with Neumann's triangle 
superimposed. 



and 0.16N/ni. There are four different regimes of be- 
haviour. In regime I, where R ^ Ts/E,h, Young's law 
holds in agreement with the analysis above. In regime II, 
as the size of the droplet shrinks towards Tg/E (though 
R ^ h), 9 starts to increase. This occurs because the 
elastic energy required to deform the substrate into a 
wetting ridge acts like a line tension, causing the droplet 
to attempt to reduce its wetted surface area. This makes 
the substrate appear less wetting. In regime IV, when 
the droplet becomes much smaller and Rq <^ "^s/E, h 
then, as we showed in Section IIIB, the substrate re- 
sponds like a fluid and the system behaves like the clas- 
sical three- fluid problem shown in Figure 1(b). Thus 9 
reduces to the value given by Neumann's triangle, with 
substrate surface tensions replaced by surface stresses. 
Unfortunately, the linear constraints of the model, de- 
tailed further in Appendix A, mean that it is difficult 
to estimate the contact angle behaviour in regime III. 
However it is interesting to note that there must be a 
pronounced maximum in 9 as R ^ 0. Note that we do 
not assume that T^ — 7^ anywhere in this calculation. 

As well as showing the change in contact angle with 
droplet size in Figure 4, we also show the free energy 
of droplet adhesion, F, changes with Rq. This is nor- 
malised by the free energy of adhesion for a rigid sub- 
strate Frigid- In regime II, there is an increase in the free 
energy above that of a rigid substrate, while in regime IV 
the free energy drops to a lower energy state than that 
in the rigid-substrate case. 

Each of the regimes above should be observable exper- 
imentally. For typical liquid droplets on gels (with kPa- 
scale moduli) , the elastocapillary length is in the micron 
scale. Thus, droplets in regime II should be readily ob- 
servable with light microscopy. Droplets in regime IV 
could be readily generated by condensation and should 
be observable with atomic force microscopy. In that case, 
our results suggest that condensation will occur more 



rapidly onto a soft substrate than onto a rigid substrate 
due to the lower energy state of droplets on the softer 
substrate, as shown in Fig. 4(b). This is in agreement 
with the experiments of Sokuler et al. [29] and the theory 
of Eslami and Elliott [12]. By contrast, the elastocapil- 
lary length for soft elastomers and rubbers {E ^ IMPa) 
is ^ 50nm, so regime I and II drops will be achievable, 
but interfacial forces will likely be important in regimes 
III and IV, which will lead to other physics being impor- 
tant in the problem. Finally the elastocapillary length 
for hard materials such as glass (GPa and above) is at, 
or below, molecular dimensions, so any droplet will be in 
regime I. 

It is worth briefly noting two points. Firstly, our re- 
sults also indicate that the predicted deviations in con- 
tact angle from Young's law also depend upon the size 
of the solid surface tension Ts, as shown in Figure 4(a). 
For larger values of Tg/ji, the changes in contact an- 
gle become small, while when Tg/^i is smaller, there can 
be quite significant changes in contact angle, especially 
when RqE/Ts ■€i 1. Physically this occurs because when 
Ts/ji is large, the liquid surface tension cannot overcome 
the solid surface tension. This means that the system ap- 
proximates that of a droplet on a rigid interface, and so 
the contact angle approaches that given by Young's law. 
On the other hand, as Tg/^i reduces, the stronger hquid 
surface tension causes surface deformations to increase, 
resulting in larger changes in 9. Secondly, we note that 
our predictions differ from those derived theoretically by 
Shanahan [30] and Leonforte and MuUer [11]. Shanahan 
[30] considered a local energy minimisation at the con- 
tact line in his derivation of the contact angle. This is in 
contrast to the global energy minimisation we perform 
to calculate 9. Leonforte and Muller primarily consid- 
ered nano droplets including intermolecular forces and 
line tension. They also performed a scaling analysis to 
approximate the competition of surface tension and elas- 
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FIG. 4. Effect of substrate deformability on (a) contact an- 
gle and (b) adsorption free energy. Change in free energy 
relative to the free energy of a droplet on a rigid substrate 
(F - F„g)/F„g. For both plots, -ji = 0.07N/m, E = SOOOPa, 
h — 20/iin. The different lines correspond to different val- 
ues of Ts; continuous line: Ts — O.OSN/m, dashed line: 
Ts = 1.2N/in, dash-dotted line: Ta — 1.6N/m. Roman nu- 
merials indicate different regimes of behaviour as described in 
the text. 



ticity and concluded that for very small droplets Young's 
law would hold, while for larger droplets Young's law 
would break down. 



B. Force balance at the contact line 

While Young's law holds for the macroscopic contact 
angle for droplets with R 3> 'yi/E,Ts/E, the micro- 
scopic behavior at the contact line is quite different. Fig- 
ure 3 shows the equilbrium of the three-phase system 
across length scales. Near the three-phase contact line 
(|r — R\ <t^ Ts/E,R), the wetting ridge forms a cusp 
with a well defined angle, as shown in Figure 3(c). In- 
terestingly, this angle is identical to the one predicted by 
Neumann's triangle with solid surface tensions T^, so lo- 
cal force balance between surface tensions is enforced at 
the tip of the ridge. An analytical proof of this result is 
given in Appendix B. This also means that the contact 
between the three phases behaves much like the contact 
between three fluids. As E —^ 0, the size of this fluid-like 
region grows, and the substrate behaviour approaches a 



completely fluid-like response. On the other hand, as E 
becomes large, the ridge height becomes small, as can be 
seen from Equation (17), and the response of the sub- 
strate returns to the flat, rigid surface of Young's law. 
Note that this picture starts to break down when the 
height of the peak is comparable to the range of interfa- 
cial forces ('^2-loA), which will effect the local behaviour 
at the contact line[27, 31]. 

Importantly, the presence of a small region where in- 
terfaces obey Neumann's triangle at the contact line en- 
sures balance in the out-of-plane component ignored by 
Young's construction. This resolves the apparent lack of 
force balance in Figure 1(a) [32, 33]. 



V. CONCLUSIONS 

We consider a sessile droplet placed on a soft sub- 
strate. We find that large droplets satisfy Young's law for 
wetting on solid substrates, while small droplets satisfy 
Neumann's construction for wetting on liquid substrates. 
The cross-over size is given by the elasto-capillary lengths 
Ts/E and 'ji/E. For droplets of all sizes, microscopic be- 
haviour near the contact line is fluid-like. At distances 
from the contact line much smaller than Ts/E, the sys- 
tem takes the form of the Neumann triangle where the 
upward force of the droplet surface tension is balanced 
by the surface tensions of the substrate. Over distances 
much greater than Ts/E the effects of the substrate sur- 
face tension diminish, and the substrate response is elas- 
tic. 

Our nrodel, which considers the case T si = Ts„, should 
be naturally extendable to treat more general surface 
stresses. To maintain local force balance at the contact 
line, we expect that the surfaces will rearrange so that 
7/1,, Tsi and T sv still obey Neumann's triangle. For large 
liquid surface tensions, the surface strain may be suffi- 
ciently large to reveal a strain-dependent contribution to 
Ts- While we have been careful to avoid equating the 
surface stresses in a solid, Tg, to the surface energy, 7^, 
equivalence of surface stress and surface energy has previ- 
ously been seen experimentally in several solid materials 
[17]. Specifically, we expect T to be isotropic and equal 
to 7 for polymer gels [10] - since their surfaces primarily 
consist of liquid solvent. 
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Appendix A: Calculating the macroscopic contact 
angle on a soft substrate 

In order to calculate the contact angle on a soft sub- 
strate for R > T s/E we use the specific case of a droplet 
on an incompressible, thin substrate with 7s; = ^sv = "is- 
Thus V ~ 1/2 and we assume i? 3> /i, as consid- 
ered in Section III A. For a hemispherical droplet in 
this limit, the height of the peak is independent of R, 
and the pressure contribution to Ff.i in Equation (20) 
vanishes as Uz{r,h) — >■ inside the droplet [10], giv- 
ing Pel = {1.5R'yf/E)f{Ts/Eh), where the function 
/ is the integral in Equation (14). When the macro- 
scopic contact angle changes from 90°, then the vertical 
component of surface tension reduces to 7; sin and so 
Eel Ri {1.5Rc-ff sin^ 9/ E)f{Ts/ Eh), where Re is the ra- 
dius of the contact area, shown in Figure 5. 




R Substrate 

c 



FIG. 5. Schematic diagram for the contact angle calculation. 

The surface energy contribution to the free energy can 
be calculated by assuming that the deflection of the sub- 
strate is small relative to the size of the droplet, so we 
have the scenario shown in Figure 5. Then, from ge- 



ometrical considerations, Fsurf — 27ri?^7/(l — cos 6*) -I- 
7ri?^(7g/ — 7s„), where Re — Rd sin. 9. Assuming that the 
total volume of the droplet is fixed so that V — 4:/3ttRq 
we find the equilibrium contact angle by numerically min- 
imising F ~ Fsurf {V, 9) -f Fei{V, 9). The results are then 
plotted in Figure 4(b) as a function of the free droplet 
radius i?o divided by Tg/E. The figure also shows the 
contact angle for small droplets in the limit Rq -C Tg/E, 
where the macroscopic contact angle is determined by 
Neumann's triangle, so that 9 = cos~^(7;/2Ts). 



Appendix B: Fluid-like behaviour at the contact line 

Here, we demonstrate that in the absence of long- 
ranged intermolecular forces, the behaviour close to the 
contact line always reverts to the Neumann triangle at 
sufficiently small length scales. We start with the equa- 
tion for the substrate profile under a droplet. Equation 
(13), and split the integral into two parts: 



Uz{r,h) 



E/T, 



ds 



ds — 



long 



short 



E/Ts 



(Bl) 



Since contributions from each wavenumber, s, separately 
satisfy the governing equations, each of the two parts 
represents a valid elastic solution. The first integral rep- 
resents the long wavelength (A > Tg/E) contribution to 
the surface profile, and is therefore smooth and cusp- free. 
The second integral represents the short wavelength con- 
tribution, and therefore contains all the details of the 
peak. For r — R '^ Tg/E, R, we find that the second 
integral reduces to 



short 



(r 



R,h) = ^ 



cos [s{r — R)] 



ds (B2) 



B/Ta 



which has a symmetric peak of width 0{Ts/E) with 
slope ±7;/2Ts either side. This is the Neumann triangle 
after linearisation for small surface gradients. Thus for 
regions of size ^ 0{Tg/E, R) the contact line appears 
the same as a three-fluid contact line. 
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